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ABSTRACT 

Context. The chemical composition of extremely metal-poor stars (EMP stars; [Fe/H]< 3) is a unique tracer of early nucleosyn- 
thesis in the Galaxy. As such stars are rare, we wish to find classes of luminous stars which can be studied at high spectral resolution. 
Aims. We aim to determine the detailed chemical composition of the two EMP stars CS 30317-056 and CS 22881-039, originally 
thought to be red horizontal-branch (RHB) stars, and compare it to earlier results for EMP stars as well as to nucleosynthesis yields 
from various supernova (SN) models. In the analysis, we discovered that our targets are in fact the two most metal-poor RR Lyrae 
stars known. 

Methods. Our detailed abundance analysis, taking into account the variability of the stars, is based on VLT/UVES spectra (R ^ 
43, 000) and ID LTE OSMARCS model atmospheres and synthetic spectra. For comparison with SN models we also estimate NLTE 
corrections for a number of elements. 

Results. We derive LTE abundances for the 16 elements O, Na, Mg, Al, Si, S, Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni, Sr and Ba, in good 
agreement with earlier values for EMP dwarf, giant and RHB stars. Li and C are not detected in either star. NLTE abundance cor- 
rections are newly calculated for O and Mg and taken from the literature for other elements. The resulting abundance pattern is best 
matched by model yields for supernova explosions with high energy and/or significant asphericity effects. 

Conclusions. Our results indicate that, except for Li and C, the surface composition of EMP RR Lyr stars is not significantly affected 
by mass loss, mixing or diffusion processes; hence, EMP RR Lyr stars should also be useful tracers of the chemical evolution of the 
early Galactic halo. The observed abundance ratios indicate that these stars were born from an ISM polluted by energetic, massive 
(25 - 40M o ) and /or aspherical supernovae, but the NLTE corrections for Sc and certain other elements do play a role in the choice of 
model. 

Key words. Stars: variables: RR Lyr - Stars: Population II - Stars: abundances - Supernovae: general - Galaxy: Halo 



1. Introduction 

The most direct information on the nature of early star forma- 
tion and nucleosynthesis in the Galaxy is pr ovided by the chem- 
ical compositi on of v ery metal - poor st ars (ICavrel et alJ I 2004: 
iFrancois et all 120071: lLai et all l2008t iBonifacio et alT l2009). 
Much of this information derives from studies of extremely 
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* Based on observations made with the ESO Very Large Telescope at 
Paranal Observatory, Chile (Large Programme "First Stars", ID 165.N- 
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metal-poor (EMP) giant stars, which are bright and display spec- 
tral lines of many elements. However, the surface chemical com- 
position of giants can be affected by products of nuclear burning 
in thei r interiors, which have b een mixed to the surface at a later 
stage (ISpiteet al.ll200ll2006l) . 

Alternative tracers exist, but have not been studied in the 
same detail as the giants. EMP turnoff stars are not suscepti- 
ble to surface mixing, but are intrinsically fainter, and the hotter 
stars may be affected by atomic diffusion, levitation and gravita- 
tional settling. In addition, many lines of heavy elements are not 
visible in spectra of EMP turnoff stars due to the higher opac- 
ity of their atmospheres. Nevertheless, the abundances observed 
in giants and dwarfs of similar metallicity generally agree well, 
although some discrepancies are seen. These are primarily as- 
crib ed to 3D and NLTE effec ts - see e.g., Bonif acio et alJ (2009) 
and lAndrievskv et al.1 d2010l) . 

Our original aim was to test whether red horizontal-branch 
(RHB) stars could be used as tracers of early Galactic nucle- 
osynthesis. To this end, we performed a detailed analysis of the 
two EMP stars CS 22881-039 and CS 30317- 056, originally dis- 
cover ed in the HK Survey of Beers et al. dBeers et al.l 1 1985L 
1 19921) and classified as RHB stars from their colours. During 
the analysis we discovered the stars to be variable in radial ve- 
locity, and subsequently that although the effective temperatures 
of both stars are clearly outs ide the fundamental red edge of 
the instability strip defined by iPreston et al.l d2006l) . both are in 
fact RR Lyrae varia bles. CS 22881-039 was studied already by 
IPreston et al.l (|2006), but from limited spectroscopic data that did 
not enable them to discover its variability. 

The paper is organised as follows: Section 2 describes our 
spectroscopic observations and Section 3 their relation to the 
pulsational phase of the stars at the times of observation. On 
this basis, Section 4 describes our derivation of stellar parame- 
ters, while Section 5 discusses the LTE abundance analysis and 
summarises our adopted NLTE abundance corrections. Section 6 
compares our abundances with earlier values for giants, dwarfs, 
and RHB stars, and Section 7 compares them with nucleosyn- 
thesis yields from recent SN models. Sections 8 and 9 present 
our discussion and conclusions. 



2. Observations and Data Reduction 

Our two targets were observed as part of our ESO Large 
Programme 'First Stars', using the VLT Kueyen telescope and 
UVES spectro graph at a resolving power R ~ 43,000, similar 
to that used by Pres ton et alJ d2006l) . The wavelength range 334- 
1000 nm is covered almost completely, using the blue and red 
cameras of UVES simultaneously. The observations and data re- 
ducti on were performed exactly as in th e rest of this programme 
(see iHill et al.ll2002t ICavrel et alJl2004l for descriptions), mak- 
ing our results directly comparable with those published previ- 
ously. 

Four spectra of CS 30317-056 were obtained in May/June 
2001. The last of these has very low S/N and was used only 
to determine the radial velocity, as was a spectrum obtained in 
August 2000 with the red camera, covering only the wavelength 
range 460 - 670nm. For CS 22881-039, only a single spectrum 
from June 2001 is available. 



3. CS 22881-039 and CS3031 7-056 as RR Lyr stars 

Following our discovery of the radial velocity variability of our 
two stars, they were also identified as photometric variables in 




phase 

Fig. 1. The radial-velocity curve of CS303 17-56, counting phase 
from maximum light. 



the database of the WASP project (A. Collier Cameron, priv. 
comm.). The periods of CS 22881-039 are P = 0.66876 and 
P = 0.74851 days, typical of RRab variables as expected for 
"red" RR Lyr stars (e.g. ISmithlll995L p. 45). The amplitude of 
CS 22881-039 is ~ 0.8 mag, typical of an RRb variable. For 
CS303 17-056 the amplitude is only ~ 0.35 mag, small for an 
RRb and more like an RRc; however, the overtone RRc pulsators 
have less asymmetric light curves with periods P < 0.5 days and 
are generally hotter than CS303 17-056. S everal RRb sta rs with 
small amplitude have been observed by dNemed 120041) in the 
metal-poor globular cluster NGC 5053, and the metallicity of 
our stars is substantially lower than that of any known globular 
cluster. 

The WASP data enable us to compute the phase <p of our 
UVES observations, counting phase from maximum light. Both 
light curves then rise steeply from phase (f> ~ 0.8 to maximum, 
with a broad minimum in the phase interval (f> ~ 0.3 - 0.8. 
Fortunately our spectra of both stars were taken near minimum 
light and thus in a relatively static phase, favourable for a spec- 
troscopic analysis ( Kolenberg et alJl2010h : the two observations 
near (f> = 0.1 and 0.9 in Fig. Q] were not used in the abundance 
analysis. 

Table [1] lists the times of the observations from which we 
derived radial velocities, and Fig. Q] shows the radial velocity 
curve of CS 30317-56. 

Only few RR Lyr stars have high-dispersion spectroscopic 
analyses, although they play an important role in the determi- 
nation of distan ces and their absolut e magnitude depends o n 
metallicity - see lLambert et al. (1996) , Sandstrom et a n j2oo"i . 
IPreston et al] d2006l) . and lPena et all (I2009I) . In fact. CS 22881- 
039 and CS 30317-056 seem to be the most metal-poor RR Lyr 
stars analysed in any detail so far. Thus, their properties are im- 
portant for our understanding of the physics of RR Lyr stars and 
the location of the instability strip in the HR diagram as a func- 
tion of metallicity. 



4. Stellar Parameters 

RR Lyr stars are pulsating variables with large variations o f 
r e ff and log g with pulsation phase: Following Pe fia et alJ (|2009), 
the change in temperature can reach 2000 K; in log g up to 1.5 
dex. It is thus very important to determine the parameters of the 
atmosphere of the stars at the exact moment of the observations. 

A rich Fe I spectrum is a rather good thermometer: The de- 
rived abundance must be independent of the excitation potential 
of the line used. In this paper we have determined T e s for each 
spectrum from the LTE excitation equilibrium of Fe I; the esti- 
mated uncertainty is about 200 K. Only the blue spectra were 
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Table 2. Atmospheric parameters for the target stars. Disregarding 15 nearly saturated Fe I lines in CS 22881-039 (22% of all Fe I 
lines) reduces the scatter in [Fe/H] by 37% (numbers in parenthesis). 



Stai- 




<x 


l°g# 


<T 


6 


cr 


[Fe/H] 


<X 


No. of lines 


Mass 


<x 




K 


K 


cgs 


cgs 


kms 1 


km s _I 




dex 


(Fe I, Fe II) 


M Q 


M G 


CS 30317-056 


6000 


200 


2.00 


0.25 


3.0 


0.1 


-2.85 


0.15 


113,4 


0.7 


0.1 


CS 22881-039 


5950 


150 


2.10 


0.3 


3.0 


0.1 


-2.75 


0.3(0.19) 


68(53), 6 


0.57 


0.12 



used in the determination, and only lines with an equivalent 
width larger than 1 pm. 

Fig-Elpresents the Fe I abundance from the three blue spectra 
of CS 303 17-056 as a function of excitation potential (Xex), com- 
puted with the same model: T eS = 6000 K, log g =2.0, [M/H]=- 
2.85 dex and microturbulent velocity, £ t = 3 kms -1 . As seen, 
for a common temperature of 6000 K there is no trend of Fe I 
abundance with the excitation potential of the lines for all the 
spectra. The star was thus rather stable during the observations, 
as expected for observations close to light minimum. 

CS 2288 1 -039 was observed only once, just before the steep 
rise in luminosity (at phase 0.86). The effective temperature de- 
rived from our spectrum is also rather low for an RR Lyr star, 
T e ff = 5950 K. The barycentric radial veloci ty is ~ 90.0 km/s. 

OSMARCS LTE m odel atmospheres dAsplund et all 1 19971: 
Gustafsson et al. 2008) and the ID LTE synthetic spectrum code 
turbospectrum (lAlvarez & Pled 1 19981) were used to derive 
st ellar parameters and abundances. The procedure is described 
in lCavrel et all d2004l) . 

The microturbulent velocity was determined by requiring 
that the derived Fe I abundance be independent of the equiva- 
lent width of the lines. 

Surface gravities (logg) were o btaine d from the ionization 
balance of Fe (as in IPreston et all (120061) 1. With r eff =6000 K 
and log g =2.0 there is good agreement between the abundances 
derived from Fe I and Fe II lines in all three spectra (Table [3j 
AS = [Fen/H] - [Fe^H] < O.ldex). 

In pulsating stars such as Cepheids or RR Lyr stars, the 
Baade-Wesselink method can also be used to determine the lu- 
minosity and thus the surface gravity of the star. However, for 
RR Lyraes it has been shown that precise infrared light c urves, 
which are n ot available for our st ars, must be used (e.g. Smith 
(119951) p. 31. rCacciari et al.l(ll992l) ). 

The final atmospheric parameters for each star are given in 
Table |2j [Fe/H] is the mean value over all measured Fe I and 
Fe II lines. 



Table 3. [Fe/H] as computed from Fe I and Fe II lines in our 
three "blue" spectra of CS 30317-56, for T eff =6000 K and 
log g =2.0. nl and n2 are the number of lines used in each de- 
termination. The phase is given for each spectrum. 



Spectrum 


[Fe,/H] 


nl 


[Fen/H] 


n2 


AE 


Phase 


151-7 


-2.89 


52 


-2.84 


9 


+0.05 


0.28 


152-4 


-2.80 


55 


-2.84 


8 


-0.04 


0.40 


152-5 


-2.82 


54 


-2.90 


7 


-0.08 


0.47 



Non-LTE effects have been observed in RR Lyr stars. If the 
luminosity is known (e.g. from a Baade-Wesselink study), the 
gravity can be derived directly. The mean abundance of Fe I is 
then found to be lower than that derived from the Fe II lines (ove- 
rionisation of iron). Following Lamb ert et al.1 (Q996), an abun- 
dance discrepancy of AE « -0.2 is seen in RR Lyr stars when 
the gravity is fixed from the luminosity; reducing of AE to 0.0 
requires a change of log g by about 0.6 dex. Thus, the true grav- 
ities (corresponding to the luminosity of the stars) are probably 
about 0.6 dex larger than the value given in Table|2] 

The masses given in Table [2] w ere estimated from th e evolu- 
tionary tr acks for [M/H ] = -2 27 bv lCassisi et alJ (120041) . As dis- 
cussed by Pres ton et al.l {2006), the changes in log g derived from 
these tracks are comparable to the observational uncertainty over 
a reasonable range in [Fe/H], and we estimate the uncertainty of 
the derived masses to b e ~ +0.1 M . 

Figure [3] shows the|Cassisi et al] (120041) evolutionary tracks 
in the log r e ff - log g plane together with our spec troscopic val- 
ues fo r CS 30317-056 and CS 22881-039. The IPreston etaD 
(2006) RHB stars are shown for comparison; as seen, CS 30317- 
056 and CS 22881-039 fall well within the range covered by 
these stars (i.e. outside the instability strip ending at the funda- 
mental red edge at log T e ff - 3.80). 
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Kiex 

Fig. 2. [Fe/H] vs. Xex for our three blue spectra of CS 30317- 
56 (152-5, 152-4 and 151-7). For clarity, the plots for 151-7 
and 152-5 have been shifted down and up by 0.8 dex, respec- 
tively. The computations were with the same OSMARCS model: 
r eff =6000 K, log g =2.0, [M/H]=-2.85 dex, £ = 3 kms" 1 . The 
mean abundance from the Fe I lines is given to the right of 
each plot. Abundances from Fe II are [Fe/H] 151-7 = -2.84 
dex, [Fe/H] 152-4 = -2.84 dex, and [Fe/H]i 52 - 5 = -2.90 dex. 
Accordingly, within the errors, the same model can be adopted 
for all spectra. 



weak (< 0.5 pm); in case of doubt, the ICavrel (Il988l) formula 
was used to assess the uncertainty of the equivalent width. 

If the Gaussian fit seemed uncertain, the equivalent width 
was remeasured by direct integration in IRAF and checked with 
an LTE synthetic spe ctrum computed with turbospectrum 
dAlvarez&Plezlll998l) . If the two derived abundances differed 
significantly, the turbospectrum result was adopted. For a few 
elements, the abundance could only be determined from the syn- 
thetic spectrum fit. 

We did not detect the Li I lines at 610.36 nm and 670.79 nm, 
even though the quality of the spectra was optimal in this region. 
The measured noise levels would allow to detect lines down to 
~ 4 mA in equivalent width, which leads to the upper abundance 
limits for Li given in Table |4] 

We have tried to measure the C abundance of CS 22881-039 
and CS 303 17-056 from a fit to the G band of the CH molecule. 
Since we found that T e ff and log g did not change between the 
three blue spectra of CS 30317-056 (section @), we coadded the 
spectra (correcting for the radial velocity variation) to improve 
the S/N ratio in the G-band region. 

The LIFBASE program of Luque & Crosley dLuqudl 19961) 
was used to compute line positions and gf-values; excita t ion en - 
ergies were taken from the line list of iJor gensen et al. (1996). 
The spectral regions free from contamination by atomic lines 
are indicated in the figure. Carbon is not detected, neither in 
CS 22881-039 nor in CS 30317-056, but the noise is large. For 
both stars the upper limit to the C abundance is log e(C) = 6.5 
(Fig.glior[C/H] < -2. 




Fig. 3. Spectroscopically determined log T e ff vs. log g and evolu- 
tionary trac ks for model masses 0.53 < M/M < 0.8 and [M/H] 
= -2.27 by ICassisi etal] d2004l) . CS 30317-056 is shown as a 
gr een triangle, C S 22881-039 as a red square, and the RHB stars 
of iPreston et all d2006) as blue x-es. 



Fig. 4. Synthetic spectrum fits to the CH band in CS 22881- 
039 and CS 30317-056. Crosses: observations; thick line: spec- 
tra computed with no carbon; thin lines: spectra computed for 
loge(C) = 6.5 and 7.2. The observed spectra are compatible with 
a complete absence of carbon; an upper limit to the abundance 
is log e(C) < 6.5. 



5. Abundance Analysis 

5. 1 . LTE abundances 

The LTE abundances were generally obtained from equivalent 
widths (listed in online TablelAl. determine d by fitting Gau ssian 
profiles to observed lines from the lists of iHill et aH (l2002 | ). For 
this, w e used the genetic algorithm f itline of [Francois et alj 
d2003l) . Lines were disregarded if they were blended or very 



We have also tried to measure the Na abundance in 
CS 22881-39 and CS 30317-56. In CS 30317-56 we observe a 
strong variation of the shape of the Na D lines with pulsation 
phase, </>: At <p ss 0.30 the D lines are relatively narrow and sym- 
metric, while they are broad and asymmetric at <p = 0.13. Only 
the first red spectrum (taken together with one of the blue spectra 
used to determine T e g , log g, and £ t ) was applied to determine 
the Na abundance. 
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In CS 22881-39 the D lines are broad and asymmetric and 
unsuitable for an abundance determination. Since we have only 
one spectrum, we could not determine whether the asymmetry is 
due to shock waves in the atmosphere just before the steep rise 
to light maximum or to a blend by interstellar D lines. 

The shape of the aluminium lines is very similar in the three 
blue spectra of CS 30317-56. Thus we have used the mean spec- 
trum to fit the line profiles. (Since one of the Al lines is in the 
wing of the Balmer line He, it is better to determine the abun- 
dance from a synthetic spectrum fit). We find [Al/Fe]=-0.77 for 
CS 22881-39 and [Al/Fe]=-0.51 forCS 30317-56. 

LTE abundances for the 16 elements O, Na, Mg, Al, Si, S, 
Ca, Sc, Ti, Cr, Mn, Fe, Co, Ni, Sr and Ba in CS 30317-056 and 
CS 2288 1 -039 are listed in Tableg]together with the NLTE abun- 
dance corrections discussed in the following section. Errors of 
abundances based on only one line have been estimated from 
the uncertainty of the equivalent width. For an easier compari- 
son with the other "First Stars" papers we ha ve adopted the same 
(ID, L TE) solar reference abundances from lGrevesse & Sauvall 
(1998), which differ by typic ally only ~ . 1 dex from the 
3D+NLTE solar abundances by Asplun d et al.ld2009l) . 



our two stars are listed in the online Table lAl Anlte varies from 
0.0 to 0.17 dex for the various Mg lines; an average value is 
-0.05 dex. 

For Silico n, Anlte is sma ll, typically from 0.0 to +0.25 dex, 
according to lShi et all (12009), who mainly consider dwarf stars; 
for the 390.5 nm line we adopt a mean value o f 0.05 dex from the 
14 stars of their Table 2. lPreston et al.l (120061) noted that silicon 
abundances for stars hotter than 5800 K decrease with increasing 
temperature and should thus not be trusted as chemical tracers. 
The temperatures of our stars are close to 5800 K, so we still 
trust our values (in LTE and NLTE), but note that a somewhat 
larger uncertainty (+ 0.1 dex) is probably in order. 

Sulphur abundances from the 921.2 nm line need a larger 
correction. With (Sh = 1), Anlte = -0.4 dex is interpolated be- 
tween the two sets of stellar para meters (T e g / l og^/IF e/HI = 
5500/2.0/-2.0 and 6500/2.0/-3.0) in lTakeda et all (120051 

For four of our nine Ca lines, Anlte has been taken from 
Mashonki na et alJ (120071) . who use a scaling factor of Sh =0.1. 
The individual corrections are listed in Table lAl and an average 
value in Table [4] 

No NLTE corrections currently exist for Ti. 



5.2. NLTE abundance correction 

NLTE effects become important when the aim is not just to com- 
pare results for different stars, but to confront the data with su- 
pernova models, which are independent of our incomplete un- 
derstanding of stellar atmospheres. 

NLTE calculations in the literature generally apply to either 
dwarfs or giants. We have therefore calculated NLTE abundance 
corrections AnlteQ f° r O and Mg for our stars and estimated 
other corrections from studies of dwarfs and giants. These cor- 
rections are discussed below and listed for all lines in online 
Table El 

Oxygen. NLTE co rrections for the O I triplet w ere computed 
using the Kiel code (ISteenbock & Holwegerl 19841) and MARCS 
mode l atmospheres. The O model atom is from lPaunzen et al.l 
(1999). The cross-sections for inelastic H I collisions were com- 
puted according to lDrawirj (11969). Setting the scaling factor for 
these cross-sections (Sh) to and 1, respectively, we obtain a 
correction of -0.1 dex in both cases. 

When we comp are these corre c tions to the more recent val- 
ues estimated from Fabbian et al. (2009) we find that our cor- 
rections are lower. iFabbian et aTl (12009) estimated NLTE correc- 
tions for RR Lyr stars 500 K warmer than our stars and obtained 
values of -0.5 and -0.35 for Sh = and 1, respectively. Their 
Fig. 8 shows how the NLTE corrections for oxygen vary with 
temperature and gravity and demonstrates that temperature has 
a much larger influence on the corrections than gravity; a differ- 
ence in I^nlte for Sh = of -0.4 dex corresponds to a difference 
of +1000 K in temperature. The warmer metal-poor stars yield 
a correction of -0.5 dex, compared to -0.1 dex for the cooler 
metal-poor turn-off stars. This would yield a correction of ap- 
proximately -0.3 dex for our stars with Sh = 0. 

a-elements. For magnesium we use precise NLTE corrections 
kindly calculated by S. And rievsky, using th e MUL TI and 
SYNTHV codes described in lAndrievskv et al.l (1201 Ol) . taking 
fine structure of the 3p 3P* levels into account, and adopting our 
oscillator strengths. The accurate results for individual lines in 

1 Anlte = log enlte - log elte- 



Odd-Z elements: Na, Al, Sc. We estimate A N lte = -0.2 and 
-0.1 dex for the Na Dl an d D2 lines in our RR Lyr stars, 
based on the calculations of (lAndrievskv et alj |2007). who ap- 
plied the same codes and methods as described for Mg above 
and Al below. For Al, we estimate Anlte = 0.7 dex for Al 
from Andrievsky et al. (2008), using a scale factor of 0.1 in the 
iDrawinl (1 19691) formula. 

No Anlte value s for Sc in metal-poor stars are available yet. 
IZhang et all {2008) suggest a correction of 0.3 dex for the Sun 
(an average for the three lines seen in our RR Lyr stars). Since 
we are dealing with Sc II, one might expect to get small NLTE 
corrections. However, due to the very different stellar parameters 
of our RR Lyr stars and the Sun and the strong dependence of the 
Sc abundance on microturbulence, we choose to not apply their 
Anlte value to our stars. 

Si-burning elements: Cr, Mn and Co. NLTE abundance correc - 
tions for Mn, and Co were ta ken f rom Bergema miet alJ (l2010h . 
Bergema nn & Gehrenl (120081) . and iBer gemann et al. (2009J), re- 
spectively. The Cr correction is estimated to be around 0.4 dex 
(M. Bergemann, priv. comm.). 

Neutron-capture elements: Sr Ba. NLTE corrections for S r 
in EMP stars are available from Mashonkin a & Gehrenl (1200 ll) . 
who find Anlte = 0.6 dex for stars of similar atmospheric 
parameters as ours, but do not list the lines used in this 
value. The correction is quite sensitive to the stellar parameter s 
(iBelvakova & Mashonkina! 1 9971: iMashonkina & Gehrenll2001l) . 
but the above value should be a good estimate fo r our s tars. 

Anlte for Ba is given by lAndrievskv et alJ d2009l) for the 
three lines at 455.4 nm (around 0.15 dex, 585.3 nm, and 649.6 
nm. Thus, only the resonance line (455.4 nm) can be corrected 
in our spectrum of CS 30317-056 and none in CS2288 1-039, so 
we prefer to retain the uncorrected LTE Ba abundances in both 
stars. 

6. Results 

In this section we compare the abundances for our two RR 
Lyr stars with the earlier results for EMP giants and dwarfs 
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Table 4. LTE abundance ratios and adopted NLTE corrections for our two RR Lyr stars, cr is the standard deviation and N the 
number of individual lines, 'ul' indicates an upper limit, '*' a value derived from a turbospectrumfit, and a ~ that the Anlte value 
is a rough estimate or only valid for some lines. 



Abundance 


CS 22881-039 




N 


Anlte 


CS 30317-056 


cr 


N 


I02 6 (Li) 


~ 0.6 


ul 


1 




- 1.0 


ul 


1 


rC/Fel 


0.69 


ul 






0.79 


ul 




TO/Fel 








-0.1 


0.9 


0.24 


1 


[Na/Fe] 


-0.13 


0.15 


2 


-0.2/-0.1 




0.13 


2 


IMa/Fel 


0.70 


0.30 


6 


-0.05 


0.48 


0.15 


6 


TAl/Fel 


-0.74 


0.05 


2 


+0.7 


-0.66 


0.15 


2 


TSi/Fel 


0.01 


0.05 


1 


-0.05 


0.05 


0.05 


1 


[S/Fe] 








-0.4 


0.78 


0.07 


1 


[Ca/Fe] 


0.29 


0.30 


5 


-0.19 


0.4 


0.13 


6 


[ScII/Fe II] 


0.11 


0.06 


3 




0.14 


0.06 


3 


[Til/Fe I] 


0.38 


0.25 


1 




0.63 


0.16 


6 


[Till/Fe II] 


0.35 


0.25 


17 




0.55 


0.12 


27 


[Cr/Fe] 


-0.32 


0.12 


4 


+0.4 


-0.18 


0.05 


5 


[Mn/Fe] 


-0.59 


0.05 


2 


+0.6 


-0.63 


0.05 


3 


[Co/Fe] 


0.14 


0.14 


2 


-+0.6 


0.27 


0.06 


2 


[Ni/Fe] 


-0.19 


0.07 


2 




-0.13 


0.05 


3 


[SrII/Fe II] 


0.00 


0.15 


1 


+0.6 


0.03 


0.10 


2 


[Ball/Fe II] 


-0.62 


0.09 


1 


-+0.15 


-0.32 


0.05 


2 



2.0 



0.5 



0.0 



-4.5 



CS 30317-056 A 



2.0 



1.0 - 



0.5 



0.0 



-4.0 



-3.5 



-3.0 



-2.5 



-2.0 



-4.5 



[Fe/H] 



i 1 1 C04* 

CS22881 -039 □ 
CS30317-056 A 
(P) CS22881 -039 * 

I 1 1 P06 X 

i 1 1 B09 • 




-4.0 



-3.5 -3.0 
[Fe/H] 



-2.5 
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Fig. 5. Left: NLTE values of [O/Fe] for our RR Lyr star s compar ed to the giants from Cavrel et"aTJ (|2004|) . Right: [Mg/ Fe], corrected 
for NLTE, compared to EMP giants dCavrel et al.l2004l) . dwarfs dBonifacio et a l2009h. as well as the RH B stars from dPreston et al.l 



2006), all have been corrected for NLTE effects by applying correctio ns from Andrievskv et al.l d2010l) . Results for CS 22881-039 
are shown as a red square (this paper) and asterisk (Preston et al. 2006); other symbols as in Fig. 2. 



Cavr el et al.1 d2004l) : [Francois et al. (200J; iBonifacio etaD 
(120091) . and RHB stars dPreston et al.| [2006). Systematic differ- 
ences in abundance scales between our results and the former 
should be minimal, because the observations and reductions 
were carried out i n exactly the same way. Larger differences 
from the results of lPreston et al.l d2006) might be expected, since 
different observational approaches, model atmospheres, codes, 
and methods were used to determine th e stellar parameters and 
abundances. Still, Prest on et al. (2006) find that differences in 
individual elemental abundances computed with MARCS or 
ATLAS model atmospheres should be less than 0.05 dex. 

The abundance ratios relative to Fe are listed in Tableland 
shown graphically in the following (see e.g. Fig. [5] [7] [8] and [9]l- 
All abundance ratios shown were obtained assuming LTE, ex- 
cept when explicitly noted as corrected for NLTE effects. Note 
that NLTE effects in Fe might offset all these abundances ra- 
tios by a small amount. Th is offset might be as large as 0.2 
dex fMashonki na et al.1 l2010) but these corrections still need in- 



clusion of many high energy levels to predict accurate abun- 
dance corrections. Generally, the figures show an abundance 
y-range of 2 dex, except for O, Sr and Ba . The two excep- 
tional stars CS 22949- 037 dCavrel et al.ll2004 and CS 29527-15 
(Bon ifacio et al.ll2009l) have been omitted in the comparisons. 
Overall, the detailed chemical composition of our two RR Lyr 
stars does not deviate significantly from the well-known mean 
trends for "normal" EMP giant and dwarf stars; we comment on 
individual element groups below. 

The element abundances derived from the three spectra 
listed in Table [3] are very similar and yield very homogeneous 
abundances (the differences in abundances from the individ- 
ual spectra are ~ 0.05 dex - per element). Hence, the abun- 
da nces are stable during the different pulsation phases (as noted 
in iKolenberg et al] d2010l) ). and should therefore be considered 
trustworthy as chemical tracers. Furthermore, the coadded spec- 
tra yield abundances agreeing within 0.05 dex with those ob- 
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tained from the single spectra, except for a few elements, where 
the difference is around 0. 1 dex. 



6.1. Lithium 

As expected, Li is not observed in our RR Lyr stars. During the 
evolution of the star along the RGB, the increasingly deep mix- 
ing has filled the external layers with the matter of deep, hot 
layers where lithium has been destroyed by proton fusion. 

6.2. The a-Elements 

As in other EMP stars, we find [a/Fe] > in our targets, with 
relatively small dispersion. We discuss the individual elements 
below. 

Oxygen 

O was only detected in CS 30317-056. Its abundance was 
determined from the O I triplet at 777. 194 nm, whi ch is affected 
by NL TE, while the O abundance for giant stars in ICavrel et al.l 
(2004) were based on the forbidden [O I] line at 630 nm, which 
is insensitive to NLTE effects, but may be affected b y 3D effects. 
Howe ver the O I 3D effects remain small for giants (ICollet et al.l 
2006). As Fig.[5]shows, especially the corrected value of [O/Fe] 
in CS 30317-056 (Sect. lBTTb a grees well w ith the mean value for 
EMP giants found by Cavre l et all (|2004|) without regard to 3D 
effects. 



the 
the 




gth [A] 



Fig. 6. The observed spectrum of CS 30317-056 around the S 
line 921.3 nm (open diamonds). The dotted blue line shows the 
spectrum of a fast-rotating star with mainly telluric lines (indi- 
cated by '*'). This confirms that the S line is indeed stellar and 
not telluric. The solid red line is an LTE synthetic spectrum fit to 
the S line, which confirms the abundance we obtain from direct 
integration. 



Mg, Si, S, Ca, and Ti 

Our Mg abundances are based on six lines in both stars, In 
CS 22881-39 the strong magnesium lines are asymmetric, which 
could indicate strong velocity fields inside the atmosphere (the 
star has been observed just before the sudden rise in luminos- 
ity). As a consequence, the Mg abundance deduced from a static 
model atmosphere cannot be completely reliable. 

Small positive NLTE corrections have been applied to all 
Mg abundances in Fig. [5] For our RR Lyr stars (see Table HJ 



y were calculated sp ecifically for this study (Sect. 15.2k for 
[Preston et all (120061) HB stars we estimate A NLTE = 0.2 - 0.3 
dex, depending on the stellar parameters. This correction may be 
overestimated, since we found a much smaller value when get- 
ting the e xact calculations compar ed to estimating the correction 
based on Andri evskv et al.l d2010l) . 

Our LTE abundances of Si are based on the single line at 
390.5 nm and exhibit large scatter, as also found by Pr eston et al.l 
(2006). Our abundances are somewhat lower than those for gi- 
ant sb^JCa^re^^l]|2003), but agree with those for turnoff stars 
bv lBonifacio et al] d2009l) . This differe nce may be lin ked to the 
temperature dependence noted by Pres ton et al] (|2006), who rec- 
ommended to only trust abundances for stars with temperatures 
below 5800 K. Our stars are close to this limit, so we expect our 
Si abundances to be fairly reliable. (The difference between the 
dwarfs and giants could be due to NLTE effects). 

Sulphur is an interesting element, because it has so far 
only been detected in a few EMP stars (e.g., CS 29497-030, 
ISivarani et all (12004 1. Here it is measurable in CS 30317-056 
only (see Fig. [6j. An LTE analysis of the strong triplet at 
921.286 nm yields [S/Fe] = 0.78 dex, abov e the mean v a lue by 
iNissen et al] d2007l) . but i n agreement withlCaffau et al.l (12005b 
and other sources listed in N issen et al] ( 2007 ). Using the value 
of Anlte = -0.4 dex from lTakeda et al] (12005 ). the agreement is 
essentially perfect. 

Our [Ca/Fe] ratios from six Ca I lines , whether der i ved in 
LTE or NLTE, a gree within the errors with lPreston et al .1 d2006h 
and lCavrel et"aTT(l2004l) . The LTE abundanc e of Ti in CS 22881 - 
039 match the Ti abundances for giants dCavrel et al] 12004). 
while tha t for CS 30317-056 matches the higher values for 
dwarfs by Bonifaci o et al. ( 20091): both a re in the range found for 
EMP RHB stars bv lPreston et al] d2006h . The Ti I and Ti II LTE 
abundances agree well, within 0.08 dex, supporting the gravity 
derived from the ionisation balance of Fe. 




-4.5 -4.0 



-3.5 -3.0 
[Fe/H] 



-2.5 



-2.0 



Fig. 7. lAl/Fel for our RR Lyr stars and the RHB stars by 
iPreston et ail C2006), as well as for giants ICavrel et all (|2004), 
all corrected for NLTE; symbols as in Fig. [5] 



6.3. The Odd-Z Elements Na, Al, and Sc 

Na and Al abundances derived from the resonance lines of the 
neutral atom s are v ery s ensitive to NLTE effec ts. According to 
ICavrel et al.l d2004h and iBonifacio et all d2009l) . A NLTE is of the 
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Fig. 8. [Cr/Fe] and [Co/Fe] vs. [Fe/H] (symbols as in Fig. [5]), 
al l derived in LTE. Th e RR Lyr stars agree well with the results 
of IPreston et all (l2006h . lCavrel et all (12004 and lBonifacio et al l 
d2009i) . 



order of -0.2/ -0. 1 dex for Na and +0.70 dex for Al. Fig.Qcom- 
pares our NLTE results for Al with those from the earlier study 
by dAndrievskv et a l. 2008). Fro m the s tellar parameters and the 
cal culations of And rievskv et al. (2008), we estimate Anlte for 
the IPreston et all (120061) RHB stars to be 0.8 - 1.0 dex for Al, 
depending on the stellar parameters. 



The values of rSc/Fe] for bo th RR Lyr stars are well within 
the range seen in ICavrel et a l. (2004). The observed trend of 
rSc/Fe] is flat, with a scatter of only 0.12 dex (Cavre Tet al.1 
2004). 



6.4. Iron-peak Elements 

To facilitate comparison with theoretical yields from supernova 
models, we discuss the two subgroups of the iron-peak elements 
separately. We note that iron itself, which is us ed as the reference 
eleme nt, may be affected by NLTE effects (Mashonkina et al.l 
l20loh . which would shift all the abundance ratios discussed in 
the following slightly up or down by a fixed amount. Most of 
the iron-peak abundance ratios are extremely tightly defined - 
so tightly, in fact, that interpretations of the slopes (notably that 
of the [Cr/Fe] relation) in terms o f metallicity-depend ent super- 
nova masses appear implausible. ICavrel et a l. (2004) discussed 
residual NLTE effects for the abundance trends as a possible 
alternative, which seems to be confirmed by the diverging abun- 



dances from Cr I and Cr II lines found bv lBonifacioetal] d~2009). 

Incomplete Silicon Burning Elements: Cr and Mn 
Our Cr abundances are based on three Cr I lines, and the Mn 
abundance on the three Mn I resonance lines near 403 nm (see 
Table [A}. Since the Mn lines are weak in these stars, no special 
treatment of the hyperfine structure was necessary. No lines of 
ionised Cr and Mn were detected in our RR Lyr stars. 

The LTE abundances of Cr and Mn for our two RR Ly r 
stars agree within 0.1 d ex with those by IPreston et alJ (12006). 
ICavrel et al.l d2004l) . and lBonifacio et all d2009t) (F ig.181 - better 
than the combined errors. iBonifacio et alJ (120091) suggest that 
NLTE effects in Cr may be significant in metal-poor stars; if so, 
[Cr/Fe] in this metallicity range should be ~ 0.1. 

Complete Silicon Burning Elements: Co andNi 
Our Co and Ni abundances are based on two neutral lines 
each in both stars. Within our estimated error, we find good 
agreement be tween our LTE abundances and those o f the 
EMP stars of ICavrel et ail d2004l) . IBonifacio etall d2009l) , and 
IPreston et alJ d2006l) ("Co only). However, we note that NLTE 
effects for Co in our RR Lyr stars are large and positive (Table|4]i. 



6.5. Neutron-Capture Elements: Sr and Ba 
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Fig. 9. [Sr/Fe] and [Ba/Fe] in EMP stars, assuming LTE (sym- 
bols as in Fig. |5J note the change in vertical scale here). Our 
RR Lyr stars agree w ell with the p revio us relations for giants , 
dwarfs and RHB st ars (ICavrel et alj2004t IBonifacio et alfeoost 
IPreston et af]|2006l) . 
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Sr and Ba are the only neutron-capture elements detected 
in our RR Lyr stars. Their abundances are determined from the 
resonance lines and hence subject to NLTE effects. However, as 
noted above, consistent NLTE corrections for Ba in our stars are 
not available in the literature, so we have chosen not to correct 
our [Ba/Fe] ratios for NLTE. 

Both [Sr/Fe] and [Ba/Fe] show huge sc atter below [Fe/H~l~ 
-3, far in excess of observational errors dPrimas et"aD t 1994: 
iMcWilliam et alj!995l:lF~rancois et alj2007l) - see Fig.|9] Taking 
this into account, our results for CS 22881-039 and CS 30317- 
056 (both corrected and uncorrected) agr ee within the com- 
bine d errors with the abu ndances found bv lPreston et alj (2006) 
and iFor & Snedenl d2010h for the RHB stars and also with the 
concordant results fo r EMP dwarfs and giants discussed by 
Bonifacio^taLl d2009l) . 



7. Supernova Models and Yields 

Having verified that our RR Lyr stars can be included in gen- 
eral samples of EMP stars, we proceed to compare the derived 
stellar abundanc es in our EMP RR Ly r stars with the supernov a 
model yields of ILimongi et al.1 (l2000t). iKobavashi et al.l (12006). 
Tominaga et alJ d2007l) and Izuta ni et al.1 (l2009h . It is generally 
believed that SNe of Type la did not contribute significantly 
to the composition of the EMP stars, so only models of core- 
collapse SNe with massive progenitor stars need to be consid- 
ered here. 

These models are characterised by a number of parame- 
ters, such as the mass of the progenitor, the explosion energy, 
the peak temperature, the mass cut, the electron density in the 
proto-neutron star, Y e 0, the amount of fallback and/or mixing, 
and the degree of anisotropy, including any pre-explosion jets. 
Variations in these parameters are manifested in characteristic 
changes in the predicted element ratios, so the observed abun- 
dance patterns can be used to constrain the values of these pa- 
rameters. For such comparisons, NLTE corrected abundances 
must of course be used whenever available - see Table [4] 

Even EMP stars are presumably the result of more than one 
SN explosion, so the assumed properties of the progenitor popu- 
lation - notably the IMF - enter implicitly into the comparison of 
models vs. observations. Only models assuming a Salpeter IMF 
have been chosen here in order to compare models on an equal 
basis without introducing differences due to changes in the IMF. 

A Salpeter IMF is the standard assumption in supernova and 
galactic chemical evolution models. Top-heavy IMF's hav e been 
tested as well by iBallero et alj d2006); Ma tteuccil d2007l) . Both 
studies find that if the presumed POP III stars more massive than 
100 M resulting in Pair Instability SNe (PISNe) with yields 
showing a strong odd/even effect, the overall impact on the early 
ISM is negligible. Only if these very massive stars kept forming 
for generations would one see an imprint of them, and this is 
not detectable in the observations. In fact, IBallero et al.l (2006) 
find that a constant or slightly varying ISM yields the best results 
when comparing to observations of metal-poor stars. 

Fig. [TO] shows that the various model yields differ most 
strongly for certain elements, which therefore become the most 
important diagnostics of the different model features. These ele- 
ments are, n otably, Na, Mg, Al, Sc, Mn and the neutron-capture 
element Sr. lLai et al.l (2008) found that Si, Ti, and Cr exhibit 
clear trends with stellar parameters such as the temperature, so 
these elements are less useful in this context. 
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2 Y e , is the so-called electron fraction, which describes the number of 
electrons per nucleon. 



Fig. 10. Observed abundance ratios in our RR Lyr stars from 
Table [4] (symbols as in Fig. |5jvs. pred icted yields fo r mod- 
els of super- an d hypernovae (HN) by: ILimongi et alJ (2000) 
(full black line ), Kobav ashi etal] (120061) (dotted: SN; dashed: 
HN), iTominaga et al.l (120071) (dot-dashed and tripl e-dot-dashed: 
20 and 40 M HN; long-dashed: 13 M Q SN) and llzutani et al.l 
(120091) 25 M HN; thick blue solid line. 



Only models that provide the fit to the abundances ob- 
served are inclu ded in Fig . [TO] Mo dels of the putative, extremel y 
massive PISNe dHeger & Woosley 2002b; Depag ne et al.ll2002l) ) 
predict a much larger odd-even effect than observed. These very 
heavy SN were therefore not considered here. 

7.1. The a -elements 

The ff-elements are generally formed in massive progenitor stars 
already before the SN explosion and thus bring information on 
the nature of the first stellar population(s). E.g., O (and C) are 
primarily produced in core helium burning in massive stars, O 
also during Ne burning, and Mg is produced in shell carbon 
burning. The paucity of iron formed by SNe II relative to the 
Solar composition causes the familiar over-abundance of the a- 
elements, typically ~0.3 dex, in metal-poor stars (see Table|4|. 

Supernova mod els iHeger & Wooslevl (2002a); 
Kobavashi et al. (2006) indicate that larger yields of even- 
Z elements, also seen as a large odd-even effect, correspond 
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to a larger (total or at least envelope) mass of the progenitor 
star. However, the magnitude of the observed odd-even effect 
depends heavily on a few relatively large NLTE abundance 
corrections, especially that for Al, so computing or estimating 
Anlte correctly is of direct importance for the conclusions that 
can be drawn on the mass of the progenitor stars. 

7.2. The Odd-Z Elements 

The odd-Z elements are generally under-abundant by ~0.2 dex 
relative to the Sun. Na and Al are produced in the progenitor 
before the SN explosion, but Sc is p roduced in explosive O an d 
Si burning, and also in Ne burning dWooslev & Weaverlll995[) . 
hence in the SN itself. 

The Sc abundance appears to be the best diagnostic of the 
explosive energy of the SN: The high (=near-Solar) [Sc/Fe], and 
to a lesser degree [Ti/Fe], ratios suggest that very energetic SNe 
(Hypernovae) were the progenitors of the EMP stars we observe 
today. We note, however, that Anlte for Sc has not yet been com- 
puted, and a large negative value would affect th is conclusion. 

Moreover, according to llzutani et al.l d2009l) . high Sc and Ti 
are tracers of a high Y e to an even larger extent than a high ex- 
plosion energy. A high-energy explosion is still required, but it 
is not the primary cause of a large production of Sc and the neu- 
tron ^^arjtur^ej^ments. Mod£kjwifhJ^tifi£ially) high values of 
Y e (llzutani et al.ll2009t iTominaga et al.l 120071) and low-density 
modifications caused by, e.g., jet-like explosions and fallback 
dTomi naga et al. 2007) fit the abundance patterns of EMP stars 
better, as indicated by the relatively large amounts of Sc and Sr. 
The caveat about NLTE effects should be kept in mind, however. 

7.3. Silicon Burning Elements 

Among the iron-peak elements, Cr and Mn are produced by in- 
complete silicon burning and Co, Ni and Zn in complete silicon 
burning during supernova explosions. Incomplete silicon burn- 
ing occurs at tempe ratures of 4 ■ 10 9 K < T pea k < 5 ■ 10 9 K 
dNomoto et alJl2005h . complete silicon burning at temperatures 
above 5 ■ 10 9 K. The observed high abundances of the complete 
Si burning elements indicate that the peak temperatures reached 
by the progenitor SNe was indeed high. At higher metallicities 
than discussed here, the production of Ni (and Fe) is, of course, 
dominated by the entirely different processes in SNe la. 

7.4. Neutron-capture Elements 

At Solar and moderately low metallicity , the abundances of S r 
and Ba are dominated by the s-process dArlandini et al.l ll999). 
In EMP stars, however, the neutron-capture el e ments should 
all be due to the r- process dSpite & Spitdll978t lTruranlll988t 
iFrancois et alJl2007h . so observations of Sr and Ba in EMP stars 
should provide useful co nstraints on the site of the r-process 
dHeger & Wooslevll2002al) . 

In past decades, core-collapse supernovae have been consid- 
ered the most promising astrophysical site for the r-proces s, but 
this scenario is facing difficulties (e.g. lWanaio et al.l |2010). and 
there is presently no consensus on the site of the r-process(es). 

Unfortunately, we only have NLTE-corrected abundances for 
Sr, and Anlte for Ba might be appreciable (see Sect. 15.2b . so our 
present ability to place strong constraints is limited. Constraints 
on Y e could be derived from abundances of light neutron-capture 
elements, i.e. elements in the range 38 < Z < 48. Slight en- 
hancements of Sc (and Ti) are consistent both with a high Y e or a 



jet-like explosion (cf. the good fit of the Tominagaet al. (2007) 
models to the data in Fig. ITOt. 



8. r eff at the RR Lyr-RHB Boundary 

One of the aims of the paper of iPreston et al.l d2006t) was to de- 
termine the temperature of the red edge of the RR Lyr instability 
strip (the "fundamental red edge": FRE in the pulsation theory). 
Horizontal branch stars cooler than the FRE do not pulsate, be- 
cause at this temperature the onset of convection disrupts the 
pulsation mechanism. This limit is difficult to predict theoreti- 
cally, however, because such a prediction requires a very good 
understanding of th e role of convection in RR Lyr pulsations. 
IPreston et alj (]2006) estimated the FRE temperature by compar- 
ing the temperature distribution of RR Lyr stars and RHB stars 
in globular clusters and in field stars. He found that in both cases 
no RR Lyr star was found below log r e fl = 3.80. Accordingly, 
he adopted for [Fe/H] < -2, logT eff (FRE) = 3.80 + 0.01 or 
T eff (FRE) = 6300+ 150 K. 

CS 22881-039 and CS 30317- 056 are both signif icantly 
cooler than the FRE as defined by IPreston et al.1 (120061) : nev- 
ertheless, they are photometric vari ables w i th lig ht curves and 
periods like those of RR Lyr stars. iNemed d2004l) found three 
RR Lyr below this limit in NGC5053 ([Fe/H] ~ -2.30), and it is 
interesting that the photometric amplitude of these st ars is small 
(< 0.6mag). In a recent paper. iFor & Snedenl (|2010) compared 
the temperature distributions of non-variable HB stars and RR 
Lyr stars and estimated the temperature of the FRE in the inter- 
val -0.8 < [Fe/H] < -2.5 be about 5900 K, in good agreement 
with our result. Hence, at the lowest metallicities the tempera- 
ture o f the FRE is probably lower than estimated by Pr eston et al.l 
(120061) : Teff(FRE) « 5900 K would be probably a better estimate. 
Moreover, the RR Lyr stars found very close to the FRE seem to 
have smaller photometric amplitudes than typical RRab. 

9. Conclusions 

Our detailed analysis of two EMP RR Lyr stars from high- 
resolution VLT/UVES spectra have yielded abundances for 16 
elements, whi ch agree very well ( generally within < 0.2 dex) 
with those by ICavrel et aTT(l2004l) . Bonifacio et all {2009), and 
Pres ton et al.l d 2006), for giants, dwarfs, and RHB stars, re: 



respec- 
tively. Due to the homogeneity of the stellar abundances from the 
individual spectra, despite the pulsation, our RR Lyr stars can be 
used as chemical tracers as well as non-variable RHB stars. 

The temperature of the fundamental red edge of the insta- 
bility strip is probab ly lower at low metallicity than obtained by 
IPreston et al.1 (|2006j); our r esults suggest a value of T^{FRE) » 
5900 K, in agreement with lFor & Snedenl (l2010l) . 

Our comparison of the abundance patterns of well- 
studied EMP stars wi t h supernova yi e lds fro m the models of 
iTominaga etail (120071) : iNomoto et all (120051) with varying Y e , 
mass cut, mixing, and fallback indicates that supernovae with 
relatively large masses (up to 40 M ) and large explosion ener- 
gies provide the best overall match to the observations. Including 
the NLTE corrections leaves much the same overall picture of the 
predecessor SNe, but lowers the estimate of the progenitor mass. 

Overall, among cu rrent models, the 40 M HN model by 
Tominaga et al.l (120071) seems to provide the best fit to our LTE 
abundance trends, even though it is not straightforward to con- 
strain the mass and energy of the previous generation of SNe 
when varying some parameters such as Y e while keeping others 
fixed. On the other hand, when we compare the NLTE corrected 
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abundances to the s ame models, the 25 M Q hypernova model by 
llzutani et al.l (|2009l) appears to give the best fit (or the 20 M e 
HN model bv lTominaga et al.l d2007l) ) , because the NLTE abun- 
dances of such elements as Na, Al, S, Sc, Ti and the iron-peak 
elements differ from the predictions of the 40 M Q HN model. 
Hence, good determination of the NLTE effects for these ele- 
ments are particularly important in comparisons with SN mod- 
els. Sc is a key element in such work, and a good NLTE analysis 
of Sc is an urgent priority. Zn is another very important diagnos- 
tic element, but is not observable in such hot EMP stars as those 
discussed here. 
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Table .1. Equivalent widths and abundances of iron lines 
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3.92 


- 


— 


402.1867 


2.76 


-0.729 


BWL 


0.61 


3.72 


- 


- 


404.5812 


1.49 


+0.280 


FMW 


8.88 


3.61 


6.84 


4.02 


406.3594 


1.56 


+0.062 


BWL 


7.92 


3.65 


5.49 


3.98 


407.1738 


1.61 


-0.022 


FMW 


7.36 


3.64 


4.80 


3.97 


413.2058 


1.61 


-0.675 


BWL 


4.88 


3.66 


1.81 


3.94 


413.4678 


2.83 


-0.649 


BWL 


0.63 


3.72 


— 


— 


414.3415 


3.05 


-0.204 


BWL 


0.75 


3.59 


— 


— 


414.3868 


1.56 


-0.511 


BWL 


5.63 


3.62 


2.59 


3.93 


418.1755 


2.83 


-0.371 


BWL 


1.43 


3.84 


- 


- 


418.7039 


2.45 


-0.548 


FMW 


1.36 


3.58 


- 


— 


418.7795 


2.42 


-0.554 


FMW 


1.56 


3.63 


- 


- 


419.1431 


2.47 


-0.666 


BWL 


1.15 


3.63 


- 


- 


419.8304 


2.40 


-0.719 


FMW 


1.63 


3.79 


— 


— 


419.9095 


3.05 


+0.155 


BWL 


1.94 


3.71 


— 


— 


420.2029 


1.49 


-0.708 


FMW 


5.68 


3.74 


— 


— 


421.0344 


2.48 


-0.928 


BWL 


0.82 


3.74 


— 


— 


421.6184 


0.00 


-3.356 


FMW 


1.80 


3.84 


- 


- 


421.9360 


3.57 


+0.000 


BWL 


0.55 


3.80 




— 


422.2213 


2.45 


-0.967 


FMW 


1.18 


3.92 


- 


- 


422.7427 


3.33 


+0.266 


BWL 


1.08 


3.60 


— 


- 


423.3603 


2.48 


-0.604 


FMW 


1.44 


3.69 


— 


— 


423.5937 


2.42 


-0.341 


FMW 


2.21 


3.60 


— 


— 


425.0119 


2.47 


-0.405 


FMW 


1.98 


3.65 


— 


— 


425.0787 


1.56 


-0.714 


BWL 


5.02 


3.66 


— 


— 


426.0474 


2.40 


+0.109 


BWL+BK 


4.18 


3.56 


— 


— 


427.1154 


2.45 


-0.349 


FMW 


2.41 


3.68 


— 


— 


427.1761 


1.49 


-0.164 


FMW 


7.87 


3.74 


— 


— 


428.2403 


2.18 


-0.779 


BWL 


1.69 


3.62 


- 


- 


429.9235 


2.42 


-0.405 


BWL 


- 


- 


0.25 


3.46 


432.5762 


1.61 


+0.006 


BWL 


7.77 


3.67 


- 


— 


435.2735 


2.22 


-1.287 


BWL 


0.83 


3.81 


— 


— 


437.5930 


0.00 


-3.031 


FMW 


3.17 


3.85 


— 


— 


438.3545 


1.49 


+0.200 


FMW 


9.00 


3.63 


— 


— 


440.4750 


1.56 


-0. 142 


FMW 


7.32 


3.63 


— 


— 


441.5123 


1.61 


-0.615 


FMW 


5.45 


3.69 


— 


— 


442.7310 


0.05 


-2.924 


BWL 


2.65 


3.67 


- 


— 


444.2339 


2.20 


-1.255 


FMW 


1.09 


3.88 


— 


— 


445.9118 


2.18 


-1.279 


FMW 


0.98 


3.82 


- 


— 


446.1653 


0.09 


-3.210 


FMW 


1.80 


3.77 


- 


- 


448.2170 


0.11 


-3.501 


FMW 


1.48 


3.97 


— 


— 


449.4563 


2.20 


-1.136 


FMW 


0.99 


3.71 


— 


— 


452.8614 


2.18 


-0.822 


FMW 


2.05 


3.75 


— 


— 


487.1318 


2.87 


-0.363 


BWL 


0.97 


3.62 


— 


— 


487.2138 


2.88 


-0.567 


BWL 


0.58 


3.60 


— 


— 


489.0755 


2.88 


-0.394 


BWL 


0.90 


3.63 


0.49 


4.13 


489.1492 


2.85 


-0.112 


BWL 


1.58 


3.60 


0.47 


3.81 


491.8994 


2.87 


-0.342 


BWL 


1.14 


3.68 


0.43 


4.01 


492.0503 


2.83 


+0.068 


BWL 


2.63 


3.69 


1.11 


4.01 


495.7299 


2.85 


-0.408 


BWL 


1.08 


3.70 


0.60 


4.21 


495.7597 


2.81 


+0.233 


BWL 


3.03 


3.58 


1.28 


3.90 


500.6119 


2.83 


-0.638 


BWL+BK 


0.62 


3.64 


- 


- 


501.2068 


0.86 


-2.642 


FMW 


1.18 


3.79 






504.1756 


1.49 


-2.203 


BWL 


0.75 


3.82 






505.1635 


0.92 


-2.795 


FMW 


0.68 


3.73 






511.0413 


0.00 


-3.760 


FMW 


0.97 


3.84 






513.9463 


2.94 


-0.509 


BWL 


0.83 


3.76 






517.1596 


1.49 


-1.793 


FMW 


1.44 


3.72 






519.1455 


3.04 


-0.551 


BWL 


0.37 


3.53 






519.2344 


3.00 


-0.421 


BWL 


0.57 


3.56 






519.4942 


1.56 


-2.090 


FMW 


0.85 


3.83 






522.7190 


1.56 


-1.228 


BWL 


3.24 


3.71 


0.77 


3.94 


523.2940 


2.94 


-0.058 


BWL 


1.60 


3.63 


0.96 


4.16 


526.6555 


3.00 


-0.386 


BWL 


0.65 


3.58 


0.28 


3.97 


526.9537 


0.86 


-1.321 


FMW 


7.42 


3.95 


2.48 


3.99 


527.0356 


1.61 


-1.339 


BWL 


2.80 


3.78 


0.51 


3.91 


532.4179 


3.21 


-0.103 


BKK 


0.66 


3.53 






532.8039 


0.92 


-1.466 


FMW 


6.08 


3.83 
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Table .1. continued. 



CS 30317-056 CS 22881-039 



A 


X 


loggf 


Ref. 


EW 


A(Fe) 


EW 


A(Fe) 


nm 


eV 






pm 


dex 


pm 


dex 


532.8532 


1.56 


-1.850 


BWL 


1.31 


3.80 


1.66 


3.96 


534.1024 


1.61 


-1.953 


BWL 


1.11 


3.87 






537.1490 


0.96 


-1.645 


FMW 


5.09 


3.84 


0.90 


3.87 


539.7128 


0.92 


-1.993 


FMW 


3.41 


3.79 






540.5775 


0.99 


-1.844 


FMW 


3.52 


3.74 






542.9697 


0.96 


-1.879 


FMW 


3.81 


3.80 


0.54 


3.86 


543.4524 


1.01 


-2.122 


FMW 


2.44 


3.80 






544.6917 


0.99 


-1.914 


BWL 


3.54 


3.82 






545.5609 


1.01 


-2.091 


BWL 


2.39 


3.76 






561.5644 


3.33 


+0.050 


BKK 


0.92 


3.65 






Fe ii 
















423.3172 


2.58 


-1.900 


av 


1.06 


3.74 






492.3927 


2.89 


-1.320 


av 


1.87 


3.74 


1.34 


4.09 


501.8440 


2.89 


-1.220 


B 


2.38 


3.78 






516.9033 


2.89 


-0.870 


FMW 


3.34 


3.64 
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Appendix A: Line Lists 
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Table A.l. The tabel includes the detected elements of CS 22881-039, compared to iPreston etal] d2006h and listed next to CS 30317-056. 
Wavelength, excitation potential (x), log gf and the obtained equaivalent width with associated abundance and uncertainty of each line is also 
given. ** means that the abundance was also found in Midas. The oscillator strengths listed for Li and Sc are values taken from VALD, but fine 
structure and hyperfine structure calculated log gf 's were actually applied to obtain the abundances listed below. A 'ul' implies that the abundances 
are only upper limits. The Mg abundances were found by line fitting applying a rotational broadning of 21 km/s - the values listed after semi-colon. 
The NLTE correction for O has been calculated by R Bonifacio, while the corrections of Mg, have been calculated by S. Andrievsky. There are two 
different values of the Mg A NLTE , first mentioned are the corrections of CS 22881-039 and last mentioned of CS 30317-056. All other estimates 
of the NLTE abundance corrections have been made based on literature (see Section l5T2l and a ~ indicates that these corrections are slightly more 
uncertain. 



Element 


A 


X 


log gf 


EW 


CS 22881-039 
Abun Eiror Abun 


Preston 
EW 


EW 


CS 30317-056 
Abun Error Abun 


Anlte 


Li 


610.37 


1.85 


0.361 


1.0 


0.6 


ul 


- 










Li 


670.79 


0.00 


-0.309 








_ 


2.2 


1 .0 


ul 




CH 


430.80 


0.16 


-1.356 




6.50 


ul 


_ 




6.50 


ul 




O 1 


777.19 


9.15 


0.369 








_ 


16.8 


6.88 


0.24 


-0.1 


Na 1 


589.0 


0.00 


0.117 




;3.46 


0.13 


_ 


67.7 


3.04 


0.06 


-0.2 


Na 1 


589.59 


0.00 


-0.184 




;3.45 


0.20 


_ 


31.4 


2.85 


0.07 


-0.1 


Mgl 


382.94 


2.71 


-0.207 


187.0 


5.88 ;5.62 


0.04 


140 


141.6 


5.06 


0.02 


0.02/0.13 


Mg 1 


383.23 


2.71 


0.146 


196.0 


5.66 ;5.62 


0.08 


_ 


171.2 


5.27 


0. 1 9 


-0.02/0.05 


Mg 1 


383.83 


2.72 


0.415 


216.0 


5.52 ;5.62 


0.03 


_ 


173.8 


4.98 


0. 10 


0.11/0.08 


Mgl 


416.37 


4.34 


-1.000 








_ 


17.70 


5.45 


0.20 




Mg 1 


435.19 


4.34 


-0.525 








_ 


28.50 


5.2 1 


0.13 




Mgl 


457.11 


0.00 


-5.393 








_ 


6.30 


5.28 


0. 10 




Mg 1 


470.30 












33 










Mg 1 


517.27 


2.71 


-0.380 


179.0 


5.56 ;5.59 


0. 1 1 


_ 








-0.01/0.03 


Mgl 


518.36 


2.72 


-0.158 


200.0 


5.72 ;5.59 


0.02 


_ 








-0.17/-0.08 


Mg 1 


552.84 


4.34 


-0.341 


23.0 


4.87 ;5.21 


0.09 


40 








0.11/0.11 


Al 1 


394.40 


0.00 


-0.640 


56.0 


2.94 


0.09 


_ 


62.2 


3.06 


0.05 


0.7 


Al 1 


396.15 


0.01 


-0.340 


86.6 


3.02 


0.07 


86 


69.4 


2.85 


0.03 


0.7 


Si 1 


390.55 


1.91 


-1.090 


122.3 


4,81,4.0** 


0.05 


1 1 1 


115.9 


4.75 


0.02 


-0.05 


S 1 


921.29 


6.53 


0.420 








_ 


51.3 


5.12 


0.07 


-0.4 


Ca 1 


422.67 


0.00 


0.240 


192.0 


4.43 


0.08 


167 


151.4 


3.76 


0.07 


0.02 


Cat 


428.30 


1.89 


-0.220 








_ 


26.2 


4. 10 


0.20 




Cat 


431.87 


1.90 


-0.210 


20.3 


3.93 


0.08 


_ 


17.4 


3.89 


0.05 




Cat 


442.54 


1.88 


-0.360 








_ 


11.4 


3.8 1 


0.05 




Cat 


443.57 


1.89 


-0.520 








_ 


12.5 


4.02 


0.07 




Cat 


445.48 


1.90 


0.260 


44.2 


3.89 


0.06 


_ 


43.1 


3.90 


0.04 




Cat 


558.88 


2.52 


0.210 








16 


18.2 


4.01 


0.06 


0.27 


Cat 


616.22 


1.90 


-0.090 


21.9 


3.78 


0.07 


_ 








0.19 


Cat 


643.91 


2.52 


0.470 


11.2 


3.47 


0.13 


_ 








0.22 


Cat 


649.38 












17 










Sc2 


424.68 


0.31 


0.240 


106.0 


0.53 


0.05 


105 


96.1 


0.37 


0.05 




Sc2 


431.41 


0.62 


-0.100 


49.9 


0.47 


0.03 


_ 


51.7 


0.46 


0. 10 




Sc2 


440.04 


0.61 


-0.536 








28 


25.3 


0.46 


0. 10 




Sc2 


441.56 


0.60 


-0.670 








_ 


22.7 


0.55 


0.20 




Sc2 


552.68 


1.77 


-0.030 








13 


13.5 


0.79 


0.09 




Ti 1 


395.82 


0.05 


-0.177 








20 


20.5 


2.94 


0. 1 8 




Ti 1 


399.86 


0.05 


-0.060 








_ 


12.5 


2.58 


0.08 




Ti 1 


453.32 


0.84 


0.476 








13 


17.6 


2.92 


0.04 




Ti 1 


498.17 


0.84 


0.50 


12.0 


2.65 


0.2 1 


_ 


11.2 


2.65 


0.08 




Ti 1 


499.11 


0.84 


0.380 








_ 


11.6 


2.78 


0.06 




Ti 1 


499.95 


0.83 


0.250 








_ 


12.4 


2.94 


0. 10 




Ti2 


375.93 


0.61 


0.270 


200.0 


3.07 


0.30 


_ 


189.6 


2.97 


0.08 




Ti2 


376.13 


0.57 


0.170 


200.8 


3.14 


0.29 


_ 


184.9 


2.96 


0.05 




Ti 2 


391.35 


1.12 


-0.410 


147.2 


3. 15 


0.20 


_ 


128.6 


2.90 


0.04 




Ti2 


401.24 


0.57 


-1.750 


53.0 


2.50 


0.13 


_ 


57.3 


2.70 


0.02 




Ti 2 


402.83 


1.89 


-0.990 


23.3 


2.53 


0.14 


_ 


28.0 


2.77 


0.03 




Ti2 


429.02 


1.16 


-0.930 


76.4 


2.5 1 


0.10 


_ 


83.7 


2.74 


0.06 




Ti2 


430.00 


1.18 


-0.490 


113.9 


2.59 


0.08 


_ 


110.1 


2.66 


0.02 




Ti 2 


430.19 












54 


51.8 


2.62 


0.04 




Ti2 


433.79 


1.08 


-0.980 


70.8 


2.41 


0.05 


_ 


59.1 


2.41 


0.05 




Ti 2 


439.41 


1.22 


-1.770 








18 


23.7 


2.79 


0.07 




Ti 2 


439.50 


1.08 


-0.510 


125.8 


2.70 


0.04 


_ 


117.6 


2.69 


0.03 




Ti2 


439.59 


1.24 


-1.970 








14 


18.7 


2.89 


0.03 




Ti 2 


439.98 


1.24 


-1.220 


46.8 


2.49 


0.04 


- 


47.7 


2.65 


0.04 




Ti 2 


441.77 


1.16 


-1.230 


41.9 


2.36 


0.03 




60.3 


2.75 


0.03 




Ti 2 


444.38 


1.08 


-0.700 


102.0 


2.51 


0.14 


101 


104.9 


2.69 


0.04 




Ti 2 


444.46 


1.12 


-2.210 










11.8 


2.79 


0.06 




Ti2 


445.05 


1.08 


-1.510 










42.3 


2.71 


0.06 




Ti 2 


446.45 


1.16 


-1.810 










26.8 


2.83 


0.04 




Ti2 


446.85 


1.13 


-0.600 


109.3 


2.57 


0.07 


106 


108.6 


2.68 


0.04 




Ti 2 


450.13 


1.12 


-0.760 


91.3 


2.46 


0.06 


99 


96.5 


2.66 


0.06 




Ti 2 


453.40 


1.24 


-0.540 


107.0 


2.57 


0.07 




107.5 


2.70 


0.04 




Ti2 


456.38 


1.22 


-0.790 


86.5 


2.52 


0.05 




85.2 


2.65 


0.02 




Ti 2 


457.20 


1.57 


-0.230 










103.6 


2.66 


0.03 




Ti2 


512.92 


1.89 


-1.300 










13.4 


2.64 


0.1 1 




Ti 2 


518.60 


1.89 


-1.370 










14.6 


2.75 


0.03 




Ti2 


518.87 


1.58 


-1.050 


34.5 


2.42 


0.10 




34.8 


2.57 


0.03 




Ti2 


522.65 


1.57 


-1.230 










25.4 


2.57 


0.04 




Ti 2 


533.68 


1.58 


-1.630 










15.2 


2.72 


0.08 




Ti 2 


541.88 












8 










V 2 


395.20 


1.48 


-0.784 








1 1 










V 2 


400.57 


1.82 


-0.522 








1 1 










Cr 1 


425.43 


0.00 


-0.110 


52.6 


2.48 


0.02 


69 


56.6 


2.59 


0.10 


0.4 
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Table A.l. continued. 











CS 22881-039 


Preston 


CS 30317-056 


Element A 


X 


log gf 


EW 


Abun Error Abun 


EW EW 


Abun Error Abun A NL te 



Cr 1 


427.48 


0.00 


-0.230 


47.6 


2.53 


0.06 


52 


55.0 


2.68 


0.10 


0.4 


Cr 1 


428.97 


0.00 


0.361 








38 


43.8 


2.66 


0.13 


0.4 


Cr 1 


520.60 


0.94 


0.020 


24.0 


2.73 


0.11 


— 










Cr 1 


520.84 


0.94 


0.160 


27.2 


2.66 


0.04 


- 










Mn 1 


403.08 


0.00 


-0.480 


37.2 


2.03 


0.16 


37 


28.8 


1.95 


0.08 


0.6 


Mn 1 




u.uu 


-u.ozu 


JU.J 


Z.Uu 


u. iy 








U.Uj 


U.D 


Mn 1 


403.45 


0.00 


-0.811 








17 


13.0 


1.88 


0.10 


0.6 


Co 1 


384.55 


0.92 


0.010 


25.9 


2.41 


0.21 




18.2 


2.28 


0.06 


-0.6 


Co 1 


399.53 


0.92 


-0.22 










12.9 


2.33 


0.06 


-0.6 


Co 1 


412.13 


0.92 


-0.320 


11.4 


2.22 


0.27 


13 


12.6 


2.40 


0.06 


-0.6 


Ni 1 


380.71 


0.42 


-1.180 


46.0 


3.26 


0.14 




42.5 


3.27 


0.07 




Ni 1 


385.83 


0.42 


-0.970 


70.3 


3.36 


0.06 




59.8 


3.29 


0.05 




Ni 1 


547.69 


1.83 


-0.890 










12.4 


3.54 


0.08 




Sr2 


407.77 


0.00 


0.170 


159.1 


0.329 


0.13 


141 


148.1 


0.17 


0.10 


0.6 


Sr2 


421.55 


0.00 


-0.145 


133.5 


0.111 


0.16 


125 


126.9 


0.02 


0.08 


0.6 


Y2 


378.87 












19 










Ba2 


455.40 


0.00 


0.17 








37 


42.4 


-1.04 


-0.14 


0.15 


Ba2 


493.41 


0.00 


-0.150 


23.8 


-1.240 


0.09 


25 











